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COORDINATE GENERATION BY CONFORMAL 
AND QUASI CONFORMAL MAPPINGS* 

by 

C. Wayne Mastin and Joe F. Thompson 
Abstract 


Many recent advances in the solution of elliptic equations have been limited 
t> regions whose boundary contours coincide with coordinate lines of the Cartesian 
i ordinate system. The reason for this lies in the fact that when an arbitrary 
curvilinear coordinate system is used, the original elliptic equation becomes 
i:;uch more complex. For example, Poisson's equation is transformed to an equation 
vith variable coefficients and a mixed derivative term. There is no added 
; nnplexity, or additional computer storage needed for thd coefficients, if an 
orthogonal coordinate system is generated from a conformal mapping. Many nuroer- 
i cal schemes are availablecfor constructing conformal mappings. Several steps 
are generally required to find the mapping betweeen an arbitrary region and a 
rectangular region. A simple finite difference scheme will be presented for 
a nstructing a conformal mapping of a rectangular region onto a simply or 
a ably -connected region. The procedure determines the module of the region, 

; 2 boundary correspondence, and the position of the interior coordinate lines. 

c coordinate generation requires the solution of a nonlinear elliptic system 
.. th oblique derivative boundary conditions, but is easily solved by any of the 
r indard iterative techniques. 

Any study of second order linear elliptic equations will introduce the 
ft t ion of quasiconformal mappings which transform the elliptic equation to 
. . conical form (i.e., the principal part of the differential operator reduces 
t c- the Laplacian). Quasiconformal mappings have been studied extensively by 
i mplex analysts, but little work has been done on the numerical construction 
quasiconformal mappings. Some preliminary theoretical and numerical results 
indicate that the above procedure for conformal mappings can be generalized to 
> instruct quasiconformal mappings. This would allow one to simultaneously fir 
t i.e boundary contours with coordinate lines and simplify the original equation 
v.en transforming to the new coordinate system. The equation, in canonical 
; . Ttr. , could possibly be solved by methods which would not be applicable to the 
* i) ir. al equation. 

Examples of coordinate systems generated from conformal and quasiconformal 
r. ppir.gs will be presented and the accuracy of the numerical scheme will be 
i scuised. For quasiconformal mappings there is need for improvement in both 
; iie r.^te of convergence of the iterative procedure and the accuracy of the 
inverted solution. 
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ABSTRACT 

Finite difference methods have been widely used in solving 
partial differential equations by numerical methods. In this paper is 
a comparison of different finite difference techniques used to solve 
Laplace's equation. Curvilinear coordinate systems are used on two- 
dimensional regions with irregular boundaries — specifically, regions 
around circles and airfoils. 

In Chapter II, truncation errors are analyzed for three different 
finite difference methods. A comparison of the false boundary method 
and two-point and three-point extrapolation schemes, used when having 
the Neumann boundary condition, is included in Chapter III, In the 
concluding chapter, the effects of spacing and nonorthogonality in 
the coordinate systems are studied. 
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TRANSFORMATION OF TWO AND THREE-DIMENSIONAL 
REGIONS 1JY ELLIPTIC SYSTEMS 


The previous status reports have emphasized progress in the analysis 
of coordinate generation methods and error in solving problems on curvi- 
linear coordinate systems. Recent emphasis has been c>n the application 
of this analysis to study the errors resulting in numerical computations. 

A compilation of numerical examples involving the solution of Laplace’s 
equation on various coordinate systems is contained in a Thesis by M. J. 
McCoy. An extension of this work to solutions of the stream function- 
vorticlty equations is in progress. In the numerical examples, the coef- 
ficients of the transformed equations have been computed both analytically 
and numerically. The analysis and numerical results were presented at a 
recent SIAM Meeting. 

An extension of the resuits on conformal mappings can also be reported, 
For years quasiconformal mappings have been of interest to mathematicians 
because they share many properties with conformal mappings and also can 
be used to reduce linear elliptic partial differentia] equations to canon- 
ical form. Applications to physical problems have been limited because 
there was no known method for constructing a quasiconformal mapping. It 
is now possible to modify the numerical scheme used for generating coor- 
dinates by elliptic systems to construct quasiconformal mappings. This 
new discovery, along with previously reported results on conformal mappings, 
will be presented at the Conference on Elliptic Problem Solvers in Santa Fe 
on June 30, 1980. The method should be of special interest to people 
working on the direct solution of elliptic problems since mixed derivative 
terms can be eliminated by quasiconformal mappings. An abstract of the 
talk is attached. The final manuscript is being prepared and will appear 
later in the proceedings of the Conference. 
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i. nrracamraioH 


Finite difference menhes have "been widely used in solving partial 
differential equations by numerical methods. In finite difference 
techniques, partial derivatives in the equation aro replaced by 
difference quotients. For example, given a differentiable function 
TJ(x,y), the first partial derivative of U with respect to x may be 
replaced by the difference of U eval\iated at two neighboring points, 
divided by the distance between these two points. This method requires 
that the region first be divided into a grid, or mesh, as shown for 
the case of a rectangle. 



The intersection of the lines in the grid are called nodes or grid 
points or mesh points. At these points we try to approximate values of 
the solution to the problem. 

For example, consider Laplace’s equation 



. 2 , 

o «j. 

— 7 - 

V 


o 


on the above rectangle whose lines are of unit width apart in each 
direction. Acsume boundary conditions arc cpecified. By replacing 
the partial derivatives in the equation by difference quotients at the 
point P, the equation is approximated by 


A.G$ IS 
JAJJE* 



♦ (p) - i( *(lf) + *(s) + *(B) + *(V) ) ** 0 

where $(P), <I>(s), 4>(E), and 4>(W) represent approximations to the 

function ♦ at P and ita neighboring pointa in the grid. 

The collection of all auch equations for all nodea gives a 

system of equations whose unknowns approximate the solution of the 

problem at the interior nodes. We may, then, set up an iterative 

scheme in order to compute an approximate value of the function at each 

» 

point. See [l] and [ 2 ] . 

We encounter a problem, however, if we UBe a rectangular grid on 
an arbitrary two-dimensional region with irregular boundaries. 



• grid points 

The grid points, or computational nodes, may not fall on the boundary 
so that we may not be able to make full use of the boundary conditions. 
When solving an equation on such a region it nay be possible to 
construct a curvilinear coordinate system such that certain coordinpte 
lineB coincide with the boundaries of the region. A curvilinear system 


1b defined to bo a finite difference grid having the property that each 
neighborhood of a node lu topologically equivalent to a rectangular 
grid In the plane , that in, the coordinate linen my be conoidered an 
level curvee of a one-to-one transformation, We can then solve our 
problem by computing a oolution to the transformed equation in the 
rectangular region. 

For an example of a curvilinear coordinate eyntem, consider the 

tranB format ion from rectangular to polar coordinates: 

x « r con P , 
y « r sin 6. 


Define the function T by (x,y) » T(r,0). If s> 0 and 0 < 0 < 2P, I is 
one-to-one except when rsfO and onto the xy-plane. 



xy-plane. 




The level lines r » constant axe mapped onto circles in the xy-plone 
and the llneb * » constant axe mapped onto epohes. 


m 






We ehall consider different typec of curvilinear coordinate systems by 
defining r and 0 as functions of f, and n • Define S as a one-to-one 
mapping of 1< f, < N, 2 1M onto 1 < r 110, 0 f, 0 1 2Iu 
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The Jacobian matrix of the transformation T o S is 

dx dx 
dr 06 

SL .21 

_or oe 

and the Jacobian = Aa long as J i 0 , the level 

lines do not intersect; and since the mapping is one-to-one in the 
region, the inverse of the transformation exists,, Further, by the 
Inverse Function Theorem, the Jacobian matrix of the inverse of the 
transformation is the inverse of the matrix of the transformation. See 

0 . 

Whe n using curvilinear systems certain difficulties may arise. 

The coordinate system can have considerable effect on the error in the 
numerical solution to the problem. Crowder and Dalton [4]) and Blottner 
and Roache [5 J have demonstrated this for the one-dimensional cbbo. E. 
de Rivas [f>] gives a study of truncation errors in the use of nonuni- 
form grids. We may readily see that coordinate line spacing rapidly 
changes or if we use an extremely nonorthogonal system, the second 
order differences of the x and y coordinate functions may become quite 
large. Therefore, when we consider the chain rule and Taylor series 
expansions, as we will do in this study, we see that our truucation 



M 

ax 

H 

3n 

M. 

& 

H 

3n 



6 
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errors can become much worae when using Buch systems, See [7J. 

In our study we shall compare several different coordinate 
systems, methods for finding values on the boundary when having the 
Neumann boundary condition, and various finite difference schemes. Lau 
[s3 has developed a finite difference method, in a manner similar to 
that to be used in thie study, for the three-dimensional case. Also 
see [S>]. 


n 

f 
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II. FINITE DIFFERENCE METHODS AND ERROR ANALYSIS 


Wo viBh to Bolve the partial differential equation 

2 2 

t ,2 9 rh . 9 A n 

v 4> = — *• + — o = 0 

9x 9y 

on a two-dimensional, doubly-connected region D, with boundary compo- 
nents and Tg, and a cut along \' and 1’^ which are coincident. Refer 
to Figure 1. We have Dirichlet boundary conditions on rg and either 

(i) Dirichlet boundary conditions on J^, or 

(ii) Neumann boundary conditions on r^. 

In order to solve this problem what we have done is to transform 

our problem from the xy-plane to the Cn -plane by a one-to-one mapping, 

# 

where our transformed plane D is a rectangular region and the boundaries 

correspond to constant coordinate lines. Since 1 ^ and are coincident 

# # # 

in D, Tj and are reentrant values in D . Refer to Figure 2. We are, 
then, able to Bolve our problem in the £n -plane. 

In this Btudy we wish to compare the accuracy of several finite 
difference methods of solving the problem, as well as to study the 
effects of nonorthogonality and rate of change of spacing of coordinate 
lines in the region. 

First we shall consider the chain rule method. Suppose a coordi- 
nate system is given in a region D of the xy-plane and f is a function 
in C^(d). Difference expressions for the first and second order partial 
derivatives of f can be obtained by transforming the region D to a 
rectangular region D and applying the chain rule. The relationship 


between the derivatives of f in the xy-plane and those in the 5n-plane 


are given below s 
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7* = 17* 77* + N ;? + 2 ^^ + Uv 7 y 2 - 

The derivatives with respect to the xy-variables may be expressed 
in terms of derivatives with respect to the f,n~varlables provided the 

Jacobian of the transformation, J ^ ~ *r^r , doeB not vanish. 

1 95 9n 9n 35 

For example, we may solve the first two equations in (l) for and . 

Higher derivatives may then be obtained by repeated applications of these 
3 f 3 f 

expressions for ~ and ~ : 

Df /ay 9f 9y 9f \/ T 

rr = hrvr -tTT r / J 
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All derivatives with respect to the f.n-variablea may be approxi- 


mated using difference operators, so we define the following expressions 
which replace the corresponding derivatives in the above equations. 


f^(P) = ( f(Q) - f(R) ) / 2 
f n (p) = ( f(S) - f(T) ) / 2 
(5) f ff (P) - f(Q) + *(H) “ 2f(P) 

f fn ( p ) = ( f(u) - f(v) + f(w) - f(x) ) / 4 

f (P) = f(s) + f(f) - 2 f(p) 

Those are merely the second order central differences on a square mesh 
of unit width in the Sn-plane. See Figure 5» The coordinates x end y 
are also functions of 6 and n, end their partial derivatives may also be 
approximated by difference expressions in a similar manner. 

From (l) and the above method of obtaining the central differences, 


then we may set up a matrix equation 
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We may solve this by 


, provided det A ^ 0. 


The determinant of A will not vanish in our case because J jt 0. 

Next we consider a Taylor series expansion of f about P. Using 

this and the difference expressions we have already developed, we get 

* 

the following matrix equation: 
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\ 


3f 

3x 
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3f 

dy 

« 

« A 

a 2 f 

3x 2 

£n 

i 

gj Gy 

H TO 

OJ M* 

nn 

_i 


3 2 f 

3y 2 


, where 


K 


+ y Ae } 

K x n y nn + y n x nn^ 

+ V + 2 ^y c ) 


H y « 

K y nn 2 
(kfs + jf ) 


x Sn y Sn * ((x V 2xx Cn ) ^rf^n ' yx €n ) K(y 

_ x nn y nn ^**nn + x n ^ ^ x nn y nn + 2x n y n^ ^nn + y n ^ 

the error terms are 0(h 5 ), h the maximum distance between P and its 

nelghborsy and 

(x 2 ) Cn ( p )=i( x (u) 2 + x(w) 2 - x(x) 2 - x(v) 2 ) 

(*y) en (P)=H x (tfMu) + *(V)y(V) - x(x)y(x) - x(v)y(v)) 
(y 2 ) Cn (*)=i(y(t0 2 + y(w) 2 - y(x) 2 - y(v) 2 ). 

We may solve by 



_-JK*S5 


When we compare the chain rule method with the Taylor series 
method, we consider the system ( 4 ) for Taylor series and the system (in) 
for the chain rule, both of which use central difference approximations 
for the partial derivatives of x and y with respect to the Sn-variableo, 
In the Taylor series expressions from ( 4 ), only the first and second 
order terms are retained, and the truncation error la O(h^) for each 
expansion. In the expressions*, for f { and f r| , if only the first order 
terms are retained, the truncation error is o(h^), These expressions, 
though, arc precisely those of the chain rule method. We see then that 
for our chain rule method, the orror terrae of f^ and f r) are O(h^), 

If we compare the last three equations in the Taylor series 
method and the chain rule method, we note that, in order to reproduce 
the Taylor series expressions from the chain rule expressions, some 
second order terms in ( 4 ) would have to he added, as well aB the error 
term which is O(h^). Bence, considering the second order terms, the 
error in these chain rule expressions is O(h^), 

Therefore, in all five equations using the chain rule, the errors 
2 

arc 0(h ). The erroro for the Taylor series method, retaining only 
first and second order terms, are all 0(h-’). It appears then that, v.'hen 
we use central difference expressions for the coefficients in these 
systems, the Taylor series expansions ohould give better results than 
the chain rule. 

The third scheme we consider uses the chain rule in the same 

manner as before, but the difference lies in our method of finding 
2 .2 a 2 

in’ "^2* H^n ' 1 — §» e-tc * I n3,tcai °f approximating these partial 
derivatives by finite difference expressions as before, we compute these 


15 


' *»* !»*«» *• 


values analytically. The accuracy of this method can also be evaluated 
from Taylor series expansions. 

For the first order central differences in the computational 
region, 


f » — 
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where only the first and second order terms are retained. 


2 »« 3 5 2 i? 


Approximations for the second order differences are obtained in a 

s 2 f . l a 4 f B , - 

« 7? + T7 7? “***" 


similar way. If f 
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a 2 f _ l / 3 4 f 
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3 4 f 


Dn 3?3n 

then, dropping all but the first and second order terms, we get the 
following estimates! 
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1/2x24 + a 4 az + a + L*4 _i ^jlz. 
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Here we observe that the difference formulae, arrived at via the chain 
rule, for the derivatives of f are accurate only if the higher order 
derivatives of the coordinate functions x and y become progressively 
□mailer* 

A comparison of the Taylor series expansion about P and the 
analytic derivatives suggests (and our study bears this out) that the 
numerical computation would be preferred due to the appearance of error 
in the first derivative terms when the coefficients are computed 
analytically. 



Ill 


ITERATION SCHEMES AND BOUNDARY TECHNIQUES 


Prom both the ohain rule and the Taylor aeriea, we hava, for 


7 4> = 0, 
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♦c 
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Then from x t y, and 41 we computo 

+C 
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.-1 
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’f,n 


nn 


to find 8 6 and 3 4 • 


8 x 


ay 


Our tranEformed equation is 

v 2 ^ = + a 2 ♦„ + a 5 +« + a 4 *r,n + a 5 ♦ 


nn 


where A 


-1 


-&J 


and a^ = . Using central differences to 


approximate V % abo re in tenaa of the -variableo, ve get 

h 4 ■ c^4>(s,r]) + Og^Cc+l^n) + Cjt(^-l,n) + c^(5,n+l) + c^(^,n-3,) 

+ Og^(?+l p n+l) + Cy|(c+lf n~l) + + c^(5-l,n-i) 

- 0. 

We will use the Successive Orcrrolaxntion iteration scheme to approxi- 
mate 4*t 

4>^U.n) » + c 2‘K^ a ' n ) + ■yU-M) + C 4 4(?,n + 1) 

+ 0^4(5, n-l) + Cg4(^+lfn+l) + c^i|>(?+l»n“l) 

+ Cg4(^~lin+l) + c^4'(C"*l»n-l)J + 

B - (r„n) 

°1 

* ^ B-l) (€,n) 

C 1 

where - -2a^ - 2a*. and « is the relaxation factor. 3o we have, then, 

* (d) * ,<- l > ♦ ■*, . 

2(a 3+ a 5 ) 

We have compared several techniques for finding values on the 
boundary. When we have the Neumann condition, - 0, on the boundary, 
we go through an iteration scheme obtained in a manner similar to that 

3 A 

of the Bcheme for the interior points. First ~ 0 ifl computed by 

the expression 


4* = v* 

3n 9 


= 0 on n = constant. 


In the first two equations of (2), if we apply the relations for 


j)f 

3X 


f\ii 

and * “ n * we obtain 

I£ = -h/ j and ^ = l!/ J 

00 that ■ (• %L f fcb ’ fe , $ ) vhora 0 = M lit ). 
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From 


i 'nco wo have — ~ 4. — ili m q 

v dx + H 9y u# 
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•i go:, the equations 
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' ion from these equatiunB and our boundary condition - i£ ii + ii it _n 

9f, ax H 3y u 

hare 


3x 


H B 21 ~H B ll) *5 + \U B 2?. “ ft M + 


3a 


-ii* L 

+ B 24 ■ It -“in) bn + (if B 25 ■ if E is) *-m ° 

r 

d l*{ + a 2 +n + v« + d 4 *f.r, + d 5 » °. 

•horr 4^ = — B^j - ~^r B^ , . The Interior boundary component will be 
~ £ with n = 1 a false boundary used for implementing the Neumann 
oun ciary condition as indicated in Figure 4 . 

So on the false boundary n * 1, in a manner similar to that used 
m obtaining a ocheme for the interior points, we have the iteration 

=h “0 * (o)( ■ + (0 ' 1) ( f . 1 ) * yrry (v f , g> 2 * n ♦ «,♦„ 
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+ d. + _ + d K d- \. 

4- /»n ?■ mj 

We have also studied th» effectn of using forward difference 
schemes for approximating values on the boundary* One method wo 
considered uses the oxpreseiona 

fff (C»n) * (-+(f, »n+ 2) + 4f ( r , ,»■+!) - >(‘ »n ))/2 

where we still have the boundary condition M- « 0. For 1% and . 

o n o x sl y 

the chain rule formula is used and the Neumann condition Mr » 0 is 

a n 

approximated by + bg <j> n * 0* Then we get tho t lire e-point extrapo- 
lation formula 

(5) ♦ (t ,n ) * a ia ia l ? - . t(Lni« , si + 2^*. ( e > /« 8 at „ « 2 . 

Another method used wao (', ,n ) a + (f, ,nH-l) - + (c,n)» Then we 
get the two-point extrapolation expression 

(6) + G;»n) “ 4 >(f,»n+l) + (r, ) / d g at t* = 2. 

Because of the difficulty in getting the partial derivatives at 
the trailing edge of the airfoil necessary in approximating values of 

there, we have used a three-point extrapolation on the trailing edge 

* 

at n =2, and i* = 1 and F» .-3 N. 

(7) +(1,2) = >(2,2) - >(3,2) + +(4,2) 

+ (N, 2) = > (H-l , 2) - > (N-2 , 2) + * (N-3 , 2 ) 

These formulas were obtained by fitting a parabola through the throe 
known values in order to approximate the points +(l,2) and j> (lf,2) on it. 




nr. COORDINATE SYSTEMS AND RESULTS 



) 


Different coordinate systems have been used in our region in order 
to compare the results when using gride with properties such as equal 
spacing; continuous changes in spacing, or an abrupt jump in spacing, 
and orthogonality or nonorthogonality. In our study we have used 

£ H o 40 and 2i n 1 M * 50 on both circles and airfoile. 

In each case, we have taken a known solution to the problem and 
compared values of this function on the grid to those of the approxi- 
mated function which we have computed. Here we provide figures 
illustrating the coordinate systems and a graph for each comparison of 
the errors of the particular methods used for the systems. The graphs 
were plotted along either K or n equal to a constant. 

The first coordinate system we will consider is an orthogonal one 
in which we have a jump in the coordinate line spacing. The lines close 
to the body are equally spaced before the jump, from radius 1 to 2.294* 
After the jump, from radius 2.294 to 10.0, though, the distance between 
consecutive coordinate lines is again uniform, this distance has been 
increased to five timeB that between the lines before the jump. This 
coordinate system is generated by the following formulas: 

We will denote r atn = constant by r(n)<- Then 
r(2) = 1 

r(n+l) » r(n) + 2.7/(M-2) for 2 in 1 M/2 
r(n+l) = r(n) + 15.5/(K-2) for M/2 l n 1 M. 

0 = 2J1 ( (5-l)/(N-l) ). 


Refer to Figures 5 and 6# In Figure 6, ye have the same system ao in 
5» except on a larger scale, and we have only shown the first 29 coordi- 
nate lines close to the body. 

In this comparison, we consider different methods for two differ- 
ent functions: 


(0 $(*»y) * x [1 + 


v) 


(ii) <j>(x,y) a h log (x 2 + y 2 ), where x * r cos 0 and y = r sin 6. 
For (i) we compare the results of chain rule and Taylor series, 
using central differences for the approximated values of the £n-partial 
derivatives of x and y, and the Neumann boundary condition on r^, 
creating a false boundary at n - 1* For a graph of these errors along 
c s 1, refer to Figure 7. The results using Taylor aeries were better 5 
than those using the chain rule. 

In (ii) we compare the chain rule and Taylor series using central 
differences for the coefficients, and the Dirichlet boundary condition. 
For a graph comparing the errors plotted along the line 5=1, reier to 
Figure 8. Again Taylor series gives better results. 

Next we will compare a nonorthogonal system with an orthogonal 
one, using the function <j>(x,y) = £ log (x 2 + y 2 ), where x = r eosO and 
y = r sin 0, In the orthogonal system, 

r = 1 + 9( (n-2)/(M-2) ) and 0 = 2 Jl( (5-l)/(N-l) ). 

Refer to Figure 9# 

For the nonorthogonal system we use 
r =5 1 + 9( (n-2)/(M-2) ) and 0max s = 511/6, Gmin = Jl/6, 


r-1 

T 


0 = Omin + (omax + omin) + 2ll ( (C-1)/(K-1) ). 


Refer to Figure 10. 
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On each of these systems we have uaed the chain rule method, 
central differences for the coefficients, and the Dirichlet boundary 


condition on r For a comparison of the error curves along' (, » 1, see 

Figure 11. The errors were constant on each circular coordinate line, 

and the results using the orthogonal syatem were better* 

For the function <f> (x,y) « x [ 1 + , we will again uaa the 

\ * +y / 

nonorthogonal system described above ( see, Figure 10 ). We will compare 
the accuracy when using a false boundary with that when using the two- 
point and three-point extrapolation formulae, (6) .and (5), respectively, 
for values on r ^ where we have the Neumann boundary condition. In each 
case we have uded the chain rule method with central difference 
expressions for the coefficients. 

At ns 1, we have taken for the x and y values, 

xfc,l) = 2x( 5 , 2 ) - x(r,3) 
yk »l) = 2y( $2) - y( c»5) 

when UBing the false boundary. For a graph along n = 2 comparing the 
errors when using these methods, refer to Figure 12, Since the error 
curves in each of these are symmetric about the line £ = 20, we have 
only shown the points for 1 <_ C £ 20. We see that both the false 

boundary and the three-point extrapolation methods give much better 

* 

results than the two-point extrapolation. 

The next coordinate system we consider iB generated bjrxer cose, 
y = r sine, 0 = 2n ( U-l)/(N-l) ), and r = 1 + 10 ( j£§J fl-tanh(^2jj. 
Refer to Figures 13 and 14. Figure 14 is a plot of the same coordinate 
system as 13 on a larger scale, so we only show the 24 lines closest 
to the body. Again we use <}>(x f y) = x ( 1 + — 

\ x'Vy/ 
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In this coordinate system we compjire the results obtained using 
the numerically approximated values of the -partial derivatives of 
x and y, or the coefficients, in both the chain rule and Taylor series 
methods with the chain rule method using the analytic derivatives of x 
and y. In each case we use the Dirlchlet boundary condition on r^. 

With both the chain rule and Taylor series where we use the 
numerical coefficients we get good results. However, when using the 

I 

analytic derivatives with the chain rule, the errors were much larger, 
liefer to Figure 15 for a graph of the errors along K — 1. 

In the previous cases, we have considered coordinate systems and 
functions on regions around clroles. Now we move to the airfoil, 
obtained by the following transformation: 


x a r cos e, y = r sin e 
xl = x(l-b), yl » y(l-b) 
x2 = xl + b, y2 = yl 



The transformation from (x2,y2) to (x,y) is the classical Joukowski 
transformation which is a conformal mapping* For this study we have 
taken b = -.1 and used the coordinate system generated by 

r = 1 + 9( (n-2)/(M-2) ) and 6 = 2JI((<- -l)/(N-l) ) 
with the function <J>(x,y) = x( 1 + — J. Refer to Figure 1 6. In each 

of the methods we will compare, we have used the three-point extrapola- 
tion formulas (?), discussed earlier, for approximating values on the 
trailing edge. 

First we compare the chain rule and the Taylor series using 
,-mnerical coefficients with the chain rule using analytically comjjuted 


25 


coefficients. In each of these, we have used the Neumann boundary 
condition, with a false boundary. For a graph of the error curves along 
the line 4 a 1, see Figure 17. For the curves graphed along 2» see 
Figure 18. Note that these error curves are symmetric about 4 = 20, so 
that in Figure 18, we have only shown these curves for 1 t; < 20. In 
both figures, we see that at most points, the Taylor series and the 
chain rule methods using approximated derivatives both give better 
results than the chain rule with analytic derivatives. 

Finally, for the same coordinate system and function $ as in the 
immediately preceding comparison, we have used the chain rule method 
with numerically approximated derivatives in comparing the false 
boundary method with the three-point extrapolation of 4 . on where ve 
have the Neumann be »dary condition. Refer to Figure 1$ for a graph 
of the error curves along n = 2. Again, since these curves are 
symmetric about 4 = 20, we have plotted them for 1 < 4 ± 20. These 

methods give equally satisfactory results. 
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y. CONCLUSION 


The objective of this study was to determine the most ac our at a 

method, of several considered, used to Bolve numerically the partial 

2 2 

differential equation v% * + ^-4- = 0. We have looked at 

3X ay 

several different coordinate systems on circles and on airfoils. In 
the comparison of the nonorthogonal system with the orthogonal one, ve 
found that the orthogonal system gave better results. 

We have also compared different schemes on the systems. When we ' 
compare Taylor series with the chain rule we see that, in general, Taylor 
series is more accurate. When considering these two methods, which use 
numerical coefficients, versus the chain rule with analytical coeffi- 
cients, we see that the errors when using the analytical, derivatives 
may be much greater than either of the first two methods. 

Also considered in this study were the false boundary, two-point 
extrapolation, and three-point extrapolation methods used when having 
the Neumann boundary condition. We found that both the falee boundary 
an:l the three-point extrapolation techniques give satisfactory results, 
each being more accurate than the two-point extrapolation method. 
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Transformed Plane 


Figure 2 
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Figure 13 
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• Chain Rule - Numerical Coefficients 


A Taylor Series - Numerical Coefficients 
B Chain Rule - Analytical Coefficients 

Figure 15 





Figure 16 


ORIGINAL PAGE is 
F POOR QUALITY 


> 40 50 

I 


— Chain Rule - Analytical Coefficients 
■ Taylor Series - Numerical Coefficients 
A Chain Rule - Numerical Coefficients 


Figure 1? 



■ Chain Rule - Analytical Coefficients 
— Taylor Series - Numerical Coefficients 
A Chain Rule - Numerical Coefficients 


Figure 10 
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